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Summary 

This study is an experimental investigation of the 
robustness of various controllers designed for the Large 
Angle Magnetic Suspension Test Fixture (LAMSTF). 
Both analytical and identified nominal models are used 
for designing controllers along with two different types 
of uncertainty models. Robustness refers to maintain- 
ing tracking performance under analytical model errors 
and dynamically induced eddy currents, while external 
disturbances are not considered. Results show that 
incorporating robustness into analytical models gives 
significantly better results. However, incorporating incor- 
rect uncertainty models may lead to poorer performance 
than not designing for robustness at all. Designing con- 
trollers based on accurate identified models gave the best 
performance. In fact, incorporating a significant level of 
robustness into an accurate nominal model resulted in 
reduced performance. This paper discusses an assortment 
of experimental results in a consistent manner using 
robust control theory. 

1 Introduction 

Robust, tracking control of the Large Angle Magnetic 
Suspension Test Fixture (LAMSTF) consists of con- 
trolling the attitude and position of a suspended rigid 
body in the presence of external disturbances and model 
uncertainties. A detailed description of the LAMSTF is 
given in [1] and the open-loop dynamic properties of the 
magnetic suspension system is given in [2]. 

The work reported herein parallels references [3]- 
[6] and is a continuation of [7]. The above references 
consider vibration attenuation and fine-pointing control 
for a stable large flexible laboratory structure. In stark 
contrast to the above passively stable flexible structure, 
the LAMSTF system is a highly unstable rigid body. 
Furthermore, the nature of the uncertainties in the two 
systems differ; the uncertainty in the system in [3]- 
[6] is mostly due to inaccurate knowledge of damping, 
frequency and mode shapes of the structural modes and 
truncated higher frequency structural modes. While for 
LAMSTF the uncertainty is mostly due to errors in the 
linearization about, the equilibrium state, an inaccurate 
knowledge of the spatial distribution of the magnetic field, 
sensor calibration errors, and errors at the plant, input, due 
to induced eddy currents. 

The underlying motivation in this study is that 
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for a. class of complex expensive systems, for example 
space vehicles, there is a. limitation on model fidelity 
due to time- varying model variations and/or limited 
amount, of data, available for sufficiently thorough system 
identification work. The novelty in this paper consists 
of determining what factors are important, in applying 
robust, control design to a. real problem. 

A series of tests was designed and conducted t.o 
investigate the effect, of Erroneous nominal models and 
erroneous uncertainty models. The effect, of erroneous 
nominal model was investigated by comparing the closed 
loop results of controllers based on analytical and iden- 
tified plants. The effect, of erroneous uncertainty model 
was investigated by comparing the closed loop results of 
nominal to a. perturbed test, which involves placing an 
aluminium ring above the suspended ma.gent. to induce a. 
particular type of eddy current.. This externally induced 
eddy current, is in addition to the existing eddy currents 
in the five iron cores of the elect.roma.gent.s which is not. 
included in the analytical model. 


2 Linear Model of LAMSTF 

The system shown in Figure 1 is the LAMSTF. It. 
consists of five electromagnets which actively suspend a. 
small cylindrical permanent, magnet.. The cylinder is a. 
rigid body and has six independent, degrees of freedom. 
Reference [7] gives a. detailed derivation of both the linear 



Figure 1: LAMSTF Configuration 


and nonlinear models. The following provides a. synopsis 
for the linear model used in the control design. 



Defining the state vector 

£ = (w 2 ,w 3 , /3 2 j3 3 , ci , v 2 , v 3 , * 2 , x 3 f (1) 

the linearized perturbed motion about the equilibrium is 
given by 

6Z = A6£ + B6p (2) 

where: the detailed analytical expressions for A and B are 
given in [7]. These expressions depend on many physical 
constants including an analytic approximation of the 
magnetic field distribution whose numerous parameters 
were obtained by numerical modeling. The variables u,:, 
fa, ty, and Xi denote the ith angular velocity of the 
cylinder with respect to the body frame; the ith Euler 
parameter relative to the inertial frame; and the ith 
translational velocity and displacement of the centroid, 
respectively. 

It can be shown that all ten states are completely 
controllable from the five coil currents. The open-loop 
system at the equilibrium state is also very unstable. For 
a detailed discussion of the physical significance of all 
modes, the interested reader is referred to [2]. 

A set of six sensors detects the five physical variables: 
pitch and yaw angles, and three displacements of the 
feentroid. Roll of the cylinder is both uncontrollable and 
unobservable. The five physical variables sensed, denoted 
as y 1 , are related approximately linearly to the state 
vector by 


signals encountered in the laboratory. On the other 
hand, model uncertainties in the loop are more difficult to 
characterize quantitatively. The model errors consideEed 
consist of two components, namely, nominal model error 
and externally induced eddy currents (due to added 
aluminium ring). No attempt was made to model the 
deviation of the nominal model from the true model 
but instead controllers based on identified and analytical 
nominal models are compared. On the other hand, errors 
due to externally induced eddy currents are heuristically 
modeled resulting in two different types of uncertainties, 
multiplicative at the input and additive. 

Figure 2 illustrates the interconnection structure. 
Due to the high speed of the real-time digital implemen- 
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y' = CS( 


(3) 


Figure 2: Interconnections in LAMSTF system 


where C is given in [7]. The measured voltage outputs 
denoted by y, are related to physical states y 1 by 

y' = [S2P] y y = [P2S] y' (4) 


tation at 800 Hz when compared to the loop bandwidth 
of approximately 50 Hz, all controllers were designed in 
the continuous domain. The small computational time 
delay and zero-order-hold effects are ignored. 


Both the sensor (physical state to output voltage) and 
actuator systems (current to forces and torques) have 
high bandwidth and are modeled as constants (not shown 
for brevity), denoted by matrices P2S and A2F. In 
summary, input to the system consists of five currents 
into five electromagnets, Sp, and the measured outputs 
are six voltage signals, y, related to the position of the 
suspended element. 

3 Controller Design 

The main objective of the controller is to stabilize 
the cylinder and track a command signal about, the 
equilibrium state. The problem is complicated by the 
ornni presence of model errors and noise. Hence, robust 
tracking is sought. This requires the specification of a 
tracking performance index and the uncertainty set for 
which the tracking performance is supposedly guaranteed. 
Of course, this guarantee is with respect to modeled 
uncertainty which itself may be uncertain. Uncertainty 
in the uncertainty model itself is a real dilemma which 
partly motivates our work. 

Controllers that meet a certain tracking performance 
are sought for a set of plants. The set. of plants is 
defined by the nominal and uncertainty models. Both 
the linearized analytical model derived from first, prin- 
ciples and an empirical model obtained from closed 
loop system identification [8] results are used as the 
nominal models. The measurement, noise and disturbance 
uncertainties are included to reflect, realistic values from 


3.1 Tracking Performance 

The tracking performance is defined by the principal gains 
such that. 

* [G e - 1 y ) -i,(/ + GA-)- 1 i;„ a ,] < 1; Vc e [0 , oo ) (5) 


where the frequency weighting function matrix, G e = 
diag(g e i, ...,g e 5 ), are given by 
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The following parameters were chosen: uj n = 2 r/s for 
pitch and yaw axes and 1 r/s for x, y, z axes, u>d = 1000 
r/s and foe = -05 for all channels. This weight, specifies 
the steady state tracking error to be within 5 % for all 
five axes up to about. 1 to 2 rad/sec and decreases by 20 
db per decade until 1000 rad/s. This frequency weighting 
results in the tracking error reaching 100%) at, a, frequency 
of about, 200 ra.d/s for both controllers. 

Note that, the above sensitivity matrix is normalized 
to unity at, the inputs and outputs by the designed 
tracking command magnitudes 

Y ma , = diag(^, .0005, .0005, .0005) (7) 

Therefore the maximum singular value corresponding to 
this scaled transfer function from command input, to 
tracking error can be interpreted as the maximum 2-norm 
error with respect to all unit, 2-norm bounded tracking 
command vectors, scaled by Y max . 









Analytical Model 

Identified Model 

-0.0000 - 0.95561 

4.0290 

-0.0000 + 0.95561 

-.3.1592 - 5.36011 

0.0000 - 7.9697i 

-3.1592 + 5 .360 li 

0.0000 + 7.96971 

-7.6822 

-9.7764 

-0.8034 -10. 29751 

9.7764 

-0.8034 +10.2975i 

57.8061 

-58.8031 - 0.42291 

-57.8061 

-58.8031 + 0.4229i 

58.7793 

60.3627 

-58.7793 

62.3473 


Table 1: Eigenvalues for analytic and identified models 


3.2 Uncertainty Models 

In this study, the main source of model errors are due 
to (1) a relatively inaccurate analytical model which is 
used as a nominal model and (2) externally induced eddy 
current effects due to the placement of an aluminum ring 
about the z axis. 

To investigate the first effect, an identified model of 
the plant without the externally induced eddy currents 
was obtained by a closed loop system identification 
technique,. Table 1 shows the poles of the analytical and 
identified models. Neither model contains transmission 
zeros. A significant difference between the analytical and 
identified model is noted, especially at lower frequencies. 
Instead of trying to generate an uncertainty model from 
these two models, separate controllers were designed 
using different nominal models. 

To investigate the second effect, two types of uncer- 
tainties were developed and tested. Only the multiplica- 
tive uncertainty description was designed to have some 
physical basis for the actual eddy current effects. In this 
way, the effect of using an incorrect uncertainty model 
can be experimentally evaluated. 

3.2.1 Additive 

The additive uncertainty model (Figure 2) maps current 
inputs to physical outputs. Assuming a nominal current 
of 1 amp input in all channels, the perturbed currents are 
scaled by the maximum anticipated physical outputs and 
multiplied by a small factor, f a = .003 and .00015, for 
significant and insignificant additive uncertainty cases. 
While the uncertainty in all other axes are assumed 
constants, the additive uncertainty in the pitch axis is 
rolled up by a first order polynomial of the form p a = 
(s + l)/(s + 100 tt ) 100tt to simulate the effect of eddy 
currents on the pitch axis. The additive uncertainty 
weight is chosen as 

W a = / a *d*a(/(p a ,l,l,l,l) (8) 

3.2.2 Multiplicative at Input 

Figure 3 shows the eddy current uncertainty modeled as 
a multiplicative type at the input. It is known that eddy 
currents are proportional to the rate of change of the 
magnetic field and therefore rate of change of actuator 
currents. This eddy current in turn produces its own 
magnetic field which interacts with the suspended magnet 
to produce perturbation forces. These perturbed forces 
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Figure 3: Input Multiplicative Uncertainty 


and moments can be modeled as a set of perturbations 
in the actuator currents. By knowing that the induced 
eddy currents affect mostly the pitch axis, the moments 
associated with the pitch axis were weighted 40 times 
more than the other axes, i.e., T = d,iag( 4, .1, .1, .1, .1). 
It is important to note that the 40 times weighting 
used is still arbitrary in the sense that it is simply our 
representation of relative dominance. This is an example 
of an uncertainty dilemma in modeling uncertainties. The 
5 by 5 unstructured uncertainty represents the unknown 
complex relationship between the time rate of change 
in the five electromagnet currents to the equivalent 
perturbed electromagnet currents. 

The pure differentiator is approximated by a proper 
transfer function uj 0 s/(s + w 0 )/. 5 X 5 where uj 0 = 1000 
rad/sec. The parameters k/Q, 0 denote the expected 
amplitude of the eddy currents, and the values k = .2 and 
.01 were chosen to represent significant and insignificant 
eddy current cases. The parameter = 100 tt rad/sec 
denotes the reference frequency with respect to which 
the magnitude of k is scaled. The constant matrices 
[,42F] and [F2,4] (not shown for brevity) denote the 
transformation from coil current to torques and forces 
on the rigid body and vice versa, 

3.2.3 Input Disturbance and Noise 

I lie input disturbances to the plant are included to 
account for bias currents and amplifier noise, a maximum 
amplitude of .1 amp is assumed in all channels over all 
frequencies. To normalize the input disturbance to a 
unit 2-norm, the diagonal weighting matrix is given by 
W w = 1/ ^5 1.5x5 • 

Although the sensing system provides a wide dy- 
namic range and is modeled as a constant matrix, [S'21 3 ] 
(not shown for brevity), there are uncertainties due to 
measurement noise, small angle/displacenrent assump- 
tions, sensor offsets from equilibrium, and calibration 
errors. In addition, the closed-loop bandwidth can be 
limited by penalizing high frequency fjontrol activity via 
assuming a fictitious large noise at high frequencies. The 
measurement noise weighting matrix chosen consists of 
W„ = G„ * [S'21 3 ] where the first term transforms assumed 
measurement noise voltages into physical displacements 
while G n = g n * Ie x 6 where 

ffn = ( S + 5000 ) * 5 °* #sd (9) 

The term N sc i = 0.04 denotes 2-sigma noise level in the 
voltage measurements observed in the laboratory. The 
output is normalized by choosing W y = 1'“^ • 

3.3 Control Designs Implemented 

Figure 4 shows the augmented plant which includes 
the performance and uncertainty weights. The p was 









Figure 4: Inputs and outputs of the augmented plant 

minimized using D-Iv iteration [9, 10, 11]. All designs 
converged in less than 5 iterations. The MATLAB 
toolbox, p.-Tools [12] was used. 

Six controllers, l\ ,, ,, p . l \ ,, ,, r . l\ lt ,n / - b ,, nir . l\ ; rn p . 
Ki mr , are designed and tested with and without the 
externally induced eddy currents, which gives a total 
of twelve experiments. The three subscripts for the 
controllers are defined as follows: the first subscript (a, *) 
denotes the analytic and identified models, respectively, 
for the nominal plant; the second subscript (m, a) de- 
notes the multiplicative and additive uncertainties; and 
the third subscript (p,r) denotes the performance and 
robustness weightings, respectively. More precisely, the 
performance weighted controller (third subscript p) is 
actually a controller with negligible built-in robustness 
while the robustness weighted controller (third subscript 
r) is designed with a significant level (compared to 
performance weighted case) of uncertainties assumed. 

The order of the controllers was reduced to ease 
real-time implementation and to remove states that are 
very weakly controllable and observable. The order of 
the controllers was reduced to approximately 30 states 
from 40 or more states. In all cases, the reduced 
controller did not significantly degrade the predicted 
robust performance. 

In the design of controllers, since; only the assumed 
level of uncertainties differ, the robust performance index, 
p, for performance controllers produced significantly 
smaller values than robustness emphasized controllers. 
All suboptimal p values for the performance controllers 
easily satisfied the constraints while the suboptimal p 
values for the robustness emphasized controllers did not 
satisfy the desired level of robust performance up to 
a factor of 2.6 for a robust controller with additive; 
uncertainties (K aar ). However, since p is defined rela- 
tive to different incorrect nominal models and incorrect 
uncertainty structures, comparing the relative p values 
attained do not necessarily translate to performance for 
the unknown real system. The identified model provides 
an accurate nominal model but the precise level and 
the uncertainty structure of the externally induced eddy 
currents are still unknown. 

4 Results 

Figure 6 shows the results from a typical tracking ex- 
periment which consists of independent step commands 
in all five axes. The step commands were 2° for pitch 
and yaw and 2 mm for x,y, and z axes. Although there 
is considerable overshoot in all axes, the damping and 
steady state errors were satisfactory. 

Since the aluminum ring predominantly affects the 
pitch axis, only pitch responses are shown in Figures 7- 
9. The simulated, nominal experiment (without added 
ring) and perturbed experiment (with ring) for controllers 


Kaap and K a ar are given in Figure 7 and similarly for 
controllers K amp and K amr in Figure 8 and controllers 
Kimp and IC imr in Figure 9. 

The simulations shown do not include the external 
source of eddy currents induced by the added aluminum 
ring so that it is only compared to the nominal experi- 
ments. Comparing performance controllers in Figures 7 
to 9 shows that the simulated response for the identified 
model (Figure 9) is much closer to the corresponding 
experiment than the analytical models (Figures 7 and 8). 
It therefore appears that the identified model is closer to 
the true nominal system than the analytical model. 

4.1 Errors in Uncertainty Model 

Controllers K aap and K aar in Figure 7 assume additive 
uncertainties while controllers K amp and K amr in Figure 
8 assume a multiplicative form. However, only the “ro- 
bust” controllers, K aar and K amr incorporate significant 
uncertainty levels. These figures show that including mul- 
tiplicative uncertainty enhances performance (cf K amr 
vs. K amp in Figure 8 for perturbed experiments) while 
including additive uncertainty degrades performance (cf 
K aar vs. K aap in Figure 7 for perturbed experiments). 
Clearly, the choice of an uncertainty model is very 
significant and it is the burden of the control designer 
to employ physical insight when applying robust control 
techniques. 

4.2 Errors in Nominal Model 

The nominal experiments for controllers Ki mp and K amp 
(Figures 8 and 9) show that controllers based on more 
accurate nominal models give significantly better perfor- 
mance in terms of overshoot and damping. However for 
robust controllers, A",; mr and K amr , the controller based 
on a more accurate nominal model gave only slightly 
better performance. This effect is amplified for the 
perturbed experiments which can be thought of as larger 
errors in both the nominal models. 

4.3 Performance and Robustness Tradeoff 

Figure 5 depicts the anticipated performance variations 
due to differences between the nominal (G 0 ) and true 
models (Gt) for differing design uncertainties (A). The 
predicted performance drops with the distance of thgtrue 
model from the nominal model. The robustness tradeoff 
involves a reduction in overall performance in order to 
guarantee reduced levels of performance for a larger set 
of plants ( A,. ) . 

All the experimental results shown in this paper can 
be explained by the generic performance and robustness 
tradeoff diagram shown in Figure 5. Three main obser- 
vations are: 

• Let G 0 be the analytic model and Go = Gt- Since 

IIApUco G 1 1 At 2 A 0 ||(x> G 1 1 A. r 1 1 oc, , A amp Could 

give poor performance while K am r should maintain 
a reduced (by tradeoff) level of performance. In 
addition, K amp is expected to be more sensitive 
(higher slope in Figure 5) to eddy current effects than 
Kamr - These predictions are all verified in Figure 8. 

• Let G 0 be the identified model and G'i = Gt- 
Since ||G'i - G 0 ||co < ||A p ||co, Kimr should give less 
performance than Ki mp . With eddy currents, there 
is only a slight reduction in performance (shaded 
region around G'i in Figure 5). Figure 9 is consistent 





Figure 5: Effect of model uncertainty on robust- 

ness /performance tradeoff 

with the above explanation. Clearly, in this case a 
performance tradeoff for robustness is unnecessary. 

• Controller K amr (Go = Gt ) should have similar 
performance as A",; mr (G'i = Gt) since both are 
robust controllers and they satisfy \\G 0 — GV||tx> < 

1 1 A r 1 1 co. This is reasonably verified by comparing 
Figures 8 and 9. Note that for K amr , the true plant, 
may lie on the shoulder (Figure 5) and is slightly 
more sensitive to eddy currents than controller A",; mr . 


5 Concluding Remarks 

Although nothing concrete has been proven, the results 
of this study provide a better understanding of the major 
factors that influence the degree of success of a robust 
control approach when applied to real problems. In 
particular, the experiments were tailored to study the 
effects of erroneous nominal and uncertainty models. The 
experimental results support the following hypothesis: 
(1) Incorporating robustness into analytical models gives 
significantly better results. However, incorporating incor- 
rect uncertainty models can lead to poorer performances 
than not designing for robustness at all. (2) Designing 
controllers based on accurate identified models gives good 
performance without robustness considerations. Finally 
it was found that the performance obtained via robust 
control techniques using only an inaccurate analytical 
model with an approximate uncertainty model is com- 
parable to the controller performance obtained by using 
an accurately identified model but without robustness 
considerations. 
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Figure 6: Tracking response for K am r and Ki mp 
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Figure 7: Experimental and simulated pitch tracking for K aap and K aar 


Simulation 


Nominal 

Experiment 


Perturbed 
Exoeri m ent 



Figure 8: Experimental and simulated pitch tracking for K amp and K amr 
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Figure 9: Experimental and simulated pitch tracking for A",; mp and A",; mr 



















